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Abstract 

We consider the effects of going beyond the approximation of a straight string 
in mesons by using a flexible flux tube model wherein a Nambu-Goto string 
bends in response to quark accelerations. The curved string is dynamically 
identical to the straight string even for ultra-relativistic mesons except for a 
small additional radial momentum. We numerically solve the curved string 
model in the case where both ends have equal mass quarks and also the case 
where one end is fixed. No approximation of non-relativistic motion is made. 
We note some small but interesting difference from the straight string. 



I. INTRODUCTION 



As the separation between quarks becomes large, their color field contracts into a string- 
like or flux tube configuration. QCD is thus thought to resemble string theory in the large 
distance regime. The Nambu-Goto-Polyakov QCD string M is an elegant model for this 
color field from which valuable insights on the nature of QCD are expected. The most 
obvious of these is the well known picture of static linear confinement. For a dynamical 
hadron the simplest assumption for the long distance color field configuration is that the 
flux tube/string always remains straight. This was the assumption of the relativistic flux 
tube model for mesons which has enjoyed phenomenological success. 

As pointed out by Nesterenko || such rigid string models cannot strictly be considered 
Nambu-Goto strings since the string cannot remain straight during dynamical motion of the 
quarks. Recently, the present authors [[| found an explicit solution to the string equations 
describing the deviation from straightness due to non-uniform end motion. The principal 
assumptions in our calculation are to treat the motion as "adiabatic" and that the deviation 
from straightness is small. The adiabatic assumption is just that the string shape depends 
on the end motion and that there is no other explicit time dependence. By direct calculation 
we find that for all physical cases the string never bends very much even for relativistic 
motion. 

An interesting conclusion of our previous work || was that almost every dynamical 
quantity for the curved string is the same as the straight string. The angular momentum, 
energy, and momentum perpendicular to the line connecting the quarks are unchanged for 
small string deformations assuming the quark distance and angular velocity are the same for 
the curved and straight strings. 

The sole difference between the dynamically curved string and the straight string is the 
component of momentum along the line connecting the quarks. From Fig. 1 the origin of 
this radial momentum is clear. It is also evident that the contribution from the outer part of 
the string will predominate due to the larger velocity of the end of the string. The net result 
will be an inward flowing radial momentum for the case shown in Fig. 1 corresponding to 
increasing angular velocity (u > 0). For a meson the total angular momentum is conserved 
and hence if the angular velocity is increasing the distance between the quarks must be 
decreasing. We conclude that the tube radial momentum due to curvature always has the 
same sign as the quark's radial momentum. In a meson with equal mass quarks the two 
halves of the string have equal but opposite radial momenta, just as the quarks do. 

Since the curved string meson is nearly identical to the rigid (straight) string meson 
similar methods are indicated for our numerical solution. The effect of the curvature induced 
radial momentum is to shift the energy level solutions downward only a few MeV. For 
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practical purposes it is quite sufficient to ignore deviation from straightness even for massless 
quarks with resulting ultra-relativistic kinematics. A principal result thus is that one can 
use the simpler straight string approximation with negligible loss in accuracy. 

The remainder of the paper is organized as follows. In Sec. [TT] we start with the string 
plus quark action and establish the equations of motion and conserved quantities. The 
solution for the curved string is reviewed in Sec. |T| and the energy, angular momentum, and 
linear momentum are computed. We show here that the curvature of the string induced by 
general quark motion does not contribute to the energy and angular momentum. The sole 
effect of the curved string is to enhance the radial momentum of the quark. The method of 
quantization of the quark-string system and the solutions scheme is covered in Sec. |IV|. Our 
numerical results and general conclusions are given in Sections M and [VI. 



II. THE STRING-QUARK ACTION 

We begin by reviewing the formalism of the Nambu-Goto-Polyakov string and spin- 
less quark actions and some expressions for conserved quantities. The action for a quark and 
anti-quark connected by a QCD string is 

S = -^J dT J da^hh ab X%X» b Vta , - ]T J dr (myf-x^). , (2.1) 

~" * * 2 = 1,2 

where h ab is the inverse of the two-dimensional metric h a b whose indices run over r and 
a, while h = det(h ab ). X^(t, a) is the string position and X^ a = d a X^. The two quark 
coordinates ) = X^ij, (jj) at the ends of the string. We use the metric signature 

conventions that h = (— , +) and r] = (— , +, +, +). 

The action ( |2.1| ) is classically equivalent to the string tension a times the worldsheet 
area swept out by the string minus the masses of the quarks times their worldline lengths. 
Variation of the Lagrange multiplier h a b determines that it is equal to the embedding metric, 
up to a multiplicative constant 

h ab = X» a X» bV ^ . (2.2) 

Variation of the string and quark positions yields the equations of motion 

(V^hh ab X$ >a = , (2.3) 
-{-) i a*/-hh° b X%={j£) %T , (2.4) 

where the canonical quark momenta are 
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<2 - 5) 

and x^ = (xf) >r . The string momentum density is 

n M (T,<7) = = -av^/i™^,^^). (2.6) 

The total energy and momentum of the string-quark system then follows 

^otal = ^ring + Quarks = H ^ ^ + £ Pf > ( 2 ' 7 ) 

17(71 i=l,2 



and the system angular momentum is 



Because of the spherical symmetry of the string and quarks system, we may consider 
motion in the (x 1 ,^ 2 ) plane. It is useful then to use helicity coordinates 

We will adopt the notation 

X = X+, (2.10) 
and therefore X* = X~. In terms of the spatial coordinate X, the above results become 

h ab = -X° a X° b + 2 Re (X;X i6 ) , (2.11) 

r<T2 

E = P° = -a daV^hh ra X° a + J2 ( m l)i ■ ( 2 - 12 ) 



The spatial momentum can be expressed in helicity form as 



V2 

-a rfav^rX.+ ^Kx),, (2.13) 

J <Tt ■ 1 r, 



where the quark coordinate Xj(r) is 

x i (r)=X(r,<7 i ). (2.14) 
Similarly, the angular momentum J 3 perpendicular to the plane of motion is 

J 3 = 2a / da^f^hh Tb \m{y,* b X) + 2 ]T (m 7 Im(x*x)) i . (2.15) 

Jen ■ i n 
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III. THE CURVED STRING WITH ONE FIXED END 



As we showed previously, it is possible to perturb about an exact straight string solution 
to obtain the shape of a string in which one end is fixed and the other has a sufficiently 
small angular acceleration. The string motion in a meson can be considered as a special case 
where the system angular momentum and energy are conserved. We begin this section by 
re-establishing this result by an abbreviated method. In || we began with the ansatz 



where F(a) is complex and F(a) and <p(t) are small. A more convenient starting point can 
be found by observing that if X(t, a = 1) = ^j=, reparametrization invariance allows us to 
set Re(-F) = 0. We can then exploit the assumption of small deviation from straightness to 
make the simpler but related ansatz 



where the overall phase is a rigid rotation plus the sum of two phases, each dependent upon 
only one of the coordinates 




(3.1) 






(3.3) 



With X M of the above form, from Eq. (|2.11j) we see that the worldsheet metric is given as 




(3.4) 



with individual components 




(3.5) 
(3.6) 
(3.7) 



h 



aRR + a 2 R 2 uf 



and volume density 





In the above expressions we have used the notation 




(3.9) 
(3.10) 
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The equations of motion for the string ( |2.3|) hold both for the temporal X° = t and 
spatial X components. In the case of X°, we have, to first order in / and 



j ± (a)aR + 'j ± (a)a 2 Ruf 



(3.11) 



Using the easily verified identity (7±(cr)i?) ji = i?((77j_(cr)) )(T , and dropping R<p terms, we 
reduce Eq. ( |3.11| ) to 



(3.12) 



Here we introduce the notation uo = 0. It can be shown [|J that the equations of motion for 
the spatial coordinate X 1 also reduce to this same equation. Integration of Eq. (|3.12|) yields 



df_ _ R 2 cb 
da v 2 . 



a 



arcsm av± 
o- 2 7_l(o-) 



1 



(3.13) 



whose solution is 



-cruj_ arcsm <rf j_ 



arcsm av j_ 



+ Cao-^x + C 2 \/l - 0"M + <tvl 



arcsm at) j_ 



(3.14) 



In the above we have defined the transverse end velocity v± = ujR. After removal of the 
rigidly rotating phase, the string coordinate can be written in terms of a radial and a trans- 
verse piece 

1 



X = V2 {XR + tXl 



which are given by 



X 



R 



R 



a . 



r> °f(p) = ^R 2 shape(a, v±) . 
R o 

Imposing the end conditions X±(a = 0) = X±(cr = 1) = 0, we find 
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(3.15) 

(3.16) 
(3.17) 



shape(a, v±) = — 



\™± (( 



arcsin t>_|_) 2 — (arcsin av±) 2 



+ a y 1 — v\ arcsin v ± — y 1 — a 2 v \ arcsin av± 



(3.18) 



6 



The string shape for v j_ <C 1 is shown in Fig. [T]. Other examples are shown in Ref. 0. 

As we pointed out earlier, one can compare the angular momentum, energy and linear 
momentum of the curved and straight string. The results are 



T 3 

u curved 

-^curved 
p 

1 _L curved 



7 3 

^straight 



aR 2 
2v~\ 



arcsm (vjj 
v± 



straight 



aR 



arcsm v± 



Pi 



straight 



aR 



v± 



1 



(3.19) 
(3.20) 
(3.21) 



Although there is no radial tube momentum for the straight string, it is non-zero for the 
curved string. This is evident from Fig. 1 and by direct calculation 



P, 



R 



aR 



uR 2 



~T v ± 

v , arcsm v ± 

± 2 



arcsm v± . 



(3.22) 



The presence of radial string momentum is the only dynamical evidence of the curved string. 

We have considered so far in this section only the heavy-light meson case where one tube 
end is fixed. Although the general straight string case was discussed previously [0 we will 
consider here the generalization to curved strings of the equal quark mass case. In Fig. 2 
we show the shape function from — 1 < a < +1 which describes the equal mass case. It is 
evident that we may immediately find the angular momentum, energy, and linear momentum 
from the preceding expressions by using the substitution 



R 

uR 



r/2 , 

v± = our/2 



(3.23) 



which leads to the following expressions for the equal quark mass case 



7 3 

u equal mass 



ar 



4v± 



arcsm v± 



v± 



-a/1 



equal mass 



ar 



arcsm v± 



P r = — ar 



v± 

2 

ujr 



1 



v± 

— arcsm v± 
2 



arcsm v± 



(3.24) 
(3.25) 
(3.26) 



Again we note that E and J 3 for the straight and curved strings are identical. In the above 
we increased E and J 3 by a factor of two because of the two string pieces. Of course we 
must add the quark contribution to each quantity to obtain the total. 
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IV. QUANTIZATION AND METHOD OF NUMERICAL SOLUTION 



We begin by reviewing the quantization of the classical expressions for the heavy-light 
straight string. From (|2.13|) , (|3.19| ) and Q2.12j ), ( p.20|) the angular momentum and energy 



are 



J A = J = m^v ± R + aR 2 F(v ± ) , (4.1) 
^ ^ arcsin v± 

E = Mn + m^ + aR , 4.2 

v± 



where Mq is the heavy quark mass and F(v±) is defined as 

1 / arcsin v± 



In the above v± = ujR and we have used the notation 



7± = / 1 = , (4-4) 



1-vl 

- = i 1 . • (4-5) 
1 - v\ - R 2 

The radial quark velocity R has to be eliminated in favor of the radial quark momentum by 
the relation 



(p^) 2 + m 2 = m 2 ( - -hi) = — — ^ (1.6) 

which leads to 



m7 = W R -/± , (4.7) 
Wr = \/(pT k ) 2 + m 2 ■ (4.8) 



From (£n|), and fljO]) we have 

J 



n W R ^ ± v ± + aRF(v ± ) } (4.9) 

JrC 

E = M Q + W Rl± + aR aTCSmV± . (4.10) 

v± 



The above equations can be quantized by making the replacements 

1 d 2 
RdR 2 



J-^+l), P 2 R ^-^R, (4.11) 



and by symmetrizing terms which depend on both R and v± with anticommutators {{A, B} = 
AB + BA). This procedure leads to the following quantum-mechanical equation for the 
angular momentum and the energy 



+ 1) 1 

-{W R , 1± v ± } + a{R,F(v ± )} , (4.12) 



R 2 

l riir , a r „ arcsinfj, ,.«„s 

H = Mq + -{W R , 7 ±} + -{R, -} ■ (4.13) 

z z v± 

As we have shown in the only difference between the straight and the curved QCD 
string is that the later has radial momentum, 

P ^ = -aR(uR 2 )g(v ± ) , (4.14) 

where we defined 



/ , r, 9~ v ± ■ \ arcsmf_L 

g{v±) = [ 1 - \Jl - vj_ - — arcsm^i. J — — A . (4.15) 



Since p R nark = Pr — PR hc , Eq. ( |4.8| ) then becomes 



Wr = V (Pfl - # c ) 2 + m2 • ( 4 -!6) 
In order to quantize the above expressions it is convenient to write 

ujR 2 = v±R — Rv± 

= l({v±,R}-{R,v±}) . (4.17) 



After promoting v± and R to quantum-mechanical operators in ( [4.17D we can use the quan- 
tum dynamical equations of motion 

v± = -i[v ±1 H], (4.18) 
R = —i[R, H] , (4.19) 

which completely specifies the quantum version of uR 2 in terms of R, v±, and Hamiltonian 
given in (|4.13|) , and allows us to quantize p^ be by choosing appropriate symmetrization 
between the non-commuting operators. This procedure gives us 

vT e = -a\{uR\\{R,g{vi-)}} ■ (4-20) 

The final step in quantizing ( 4.16|) is to expand (pr — p^ be ) by writing 
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(PR - pT e ? =Pl- {PR, PR hc } + (PR hc ? • (4-21) 
Here, p R = and p^ he is denned by ( P7|) through ( [OUfl . 

The solution method for the straight string equations is discussed in detail in Ref. [|J. The 
fundamental problem is that the unknown operator v± cannot be analytically eliminated, so 
one must eliminate it numerically. This can be accomplished by expanding the radial wave 
function in terms of a complete set of basis states, which reduces ( }4.12| ) to a transcendental 
infinite dimensional matrix equation. By truncating the number of the basis states to a 
finite number N, the matrix equation becomes finite dimensional, and can be solved by 
iteration procedure described in [fj]. Once the matrix for v± is known, one can calculate 
the Hamiltonian matrix from (|4.13| ) and solve for the energy eigenvalues and corresponding 
wave functions. The modified algorithm for solving the curved string equations involves 
continually updating the Wr matrix during the iteration procedure in which the matrix 
v± is found. 

Note that truncation introduces the dependence of the eigenvalues on the variational 
parameter f3 describing the basis states. However, if a sufficiently large number of basis states 
are used, this dependence will exhibit a "plateau" indicating a stable solution. Examples of 
this will be shown in Sec. |V[ 

The equal quark mass meson can be quantized by a simple modification of the above. In 
this case we start with Eqs. (|3.24|) - (|3.26|) to obtain expressions analogous to ( f4.12|) , (|4.13| ), 
and O) 



£±H = {W r ,y±v±} + l{r,F{v x )} , (4.22) 

f" i a r arcsint>_i_ ^ 

H = {Wr,j±} + 2l r ' } , (4-23) 

Pt ubc = -^{^ 2 ^{^^K)}} • (4-24) 

The numerical solution is step by step the same as for the heavy-light case. 

The more general case of unequal arbitrary masses at the ends is one step more involved. 
This case was solved in || for the straight string. For mesons with unequal mass quarks, 
one must require that P| ot = 0. This is equivalent to solving for the center of momentum 
in addition to the angular momentum and energy equations. The curved string dynamical 
solution should proceed similarly to the two special cases discussed previously. 



V. NUMERICAL RESULTS 



For zero angular momentum states we observe from fl4.12j ) and fl4.22[) that v± = as 
expected. The tube radial momentum is zero and the Hamiltonians ( [4.13|) and fl4.23j ) become 
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equivalent to "spinless Salpeter" equations, which have well defined numerical solutions [0. 
We will therefore confine ourselves to non-zero angular momenta. 

The results for equal mass and heavy-light mesons are qualitatively similar and we will 
discuss these string configurations in parallel, pointing out the similarities and differences. 
As mentioned in the preceding section the set of basis states have a scale parameter (3. For 
a given truncation to iV basis states it is customary to plot energy eigenvalues as a function 
of P to identify a "plateau" indicating a stable solution. In Fig. |] (heavy-light) and Fig. [| 
(equal mass) we show the j3 dependence of the excitation energy of the p-wave ground state 
meson with massless light quarks. The energy calculation is done twice; first with the straight 
string model and then with the curved string. We note that as the number of basis states 
increases definite plateaus develop. 

The second observation we make from Figs. and f| is that the curved string energy is 
slightly less than the straight string energy. We believe this results from the reduction of the 
effective string tension due to the transverse motion of the string during the quark's motion. 
The difference in the p-wave case, shown in Fig. ^|, is just a few MeV but this varies with 
angular momentum and radial excitation. 

In Figs. H and |6] we illustrate how the curved and straight string differ for a variety of low 
lying states. The cases shown are for zero mass light quark(s) at the end(s) of the string. 
The general trend is that the difference increases from zero at i = (which must be the 
case) to about 15 MeV at the higher angular momentum values shown. 

The dependence of the energy difference for p-wave solutions on the mass of the quark 
at either end of the string is illustrated in Fig. [7| (heavy-light) and Fig. || (equal mass). The 
energy difference AE is seen to decrease rapidly as the quark mass increases. As expected, 
the curvature effect on the solution is largest for massless quarks. 

Finally, in Fig. |9] we show the Regge pattern for the curved string with one massless 
quark and the other end fixed. We exhibit the corresponding Regge pattern for massless 
quark ends in Fig. On this scale they are similar to those of the straight string. One 
observes towers of nearly degenerate states of even (or odd) angular momentum. The two 
cases are analogous except that the heavy-light Regge slope is double the light-light slope. 
Except for the value of the slope, this pattern is very similar to that of scalar potential 
confinement ||[7|- 

VI. CONCLUSION 

In this paper we have investigated the dynamical effects of the small radial momentum 
developed by a QCD string as it bends in response to the angular acceleration of the quarks 
on the ends. Our principal results are that one can successfully implement the curved string 
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solutions into a dynamical meson model with massive quarks at the string ends. The solution 
is self-consistent in that for a physical state the deviation from straightness is small. For such 
small transverse string motion the energy, angular momentum, and transverse momentum 
are identical to those of a straight string with the same angular velocity. When the small 
curvature induced radial momentum is incorporated into the numerical method one finds 
solutions that are very similar to those of the straight string. 

The curvature effects are largest for massless quarks. Even in this case one finds shifts of 
at most a few MeV. Our results indicate that for most purposes the simpler straight string 
model gives nearly correct results and if desired, more accuracy can be readily achieved. 
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FIG. 1. The solution to the string shape ( 3.17| ) with positive Co and v± <C 1. The transverse and 
radial momenta of two elements are depicted. The contribution at the larger radius will predominate 
giving a net radial momentum parallel to the quark's radial momentum. 




v, = 1 



FIG. 2. The solution to the string shape ( 3.17|) with positive Co for a string with equal mass 
quarks at its ends. The "shape function" fl3.18| ) is shown for non-relativistic and ultra-relativistic 
quark rotation. 
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FIG. 3. Dependence on the variational parameter (3 of the ground state p-wave excitation energy 
of the heavy-light Hamiltonian with the straight string (upper curves) and with the curved string 
(lower curves). We used a = 0.2 GeV 2 and light quark mass m = 0. Results shown were obtained 
with 10, 20 and 30 basis states. 
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FIG. 4. Dependence on the variational parameter (3 of the ground state p-wave excitation energy 
of the equal mass Hamiltonian with the straight string (upper curves) and with the curved string 
(lower curves). We used a = 0.2 GeV 2 and quark mass m = 0. Results shown were obtained with 
10, 20 and 30 basis states. 
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/ 

FIG. 5. Difference between the energies of the straight and curved heavy-light Hamiltonian as 
a function of the orbital angular momentum. We used a = 0.2 GeV 2 and light quark mass m = 0. 
Results shown are for the ground state (n = 1) and for the first two radially excited states. 
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FIG. 6. Difference between the energies of the straight and curved equal mass Hamiltonian as 
a function of the orbital angular momentum. We used a = 0.2 GeV 2 and light quark mass m = 0. 
Results shown are for the ground state (n = 1) and for the first two radially excited states. 
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FIG. 7. Difference between the p-wave it = 1) energies of the straight and curved heavy-light 
Hamiltonian as a function of the light quark mass. We used a = 0.2 GeV 2 , and results shown are 
for the ground state (n = 1) and for the first four radially excited states. 
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FIG. 8. Difference between the p-wave (£ = 1) energies of the straight and curved equal-mass 
Hamiltonian as a function of the quark mass. We used a = 0.2 GeV 2 , and results shown are for 
the ground state (n = 1) and for the first four radially excited states. 
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FIG. 9. Regge trajectories for the heavy-light Hamiltonian with the curved string having a 
massless quark. The vertical line illustrates the tower structure of nearly degenerate states of 
alternating angular momentum. 
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FIG. 10. Regge trajectories for tphe equal-mass Hamiltonian with the curved string having 
massless quarks at its ends. The vertical line illustrates the tower structure of nearly degenerate 
states of alternating angular momentum. 
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